NUCLEAR SCIENCE AND TECHNIQUES 34, 203 (2023) 


Simulation study of BESIII with stitched CMOS pixel detector using ACTS* 


Yi Liu,! Xiao-Cong Ai," ! Guang-Yan Xiao,” Ya-Xuan Li? Ling-Hui Wu,* Liang-Liang Wang,* Jia-Ning 
Dong,’ Ming-Yi Dong, ^ Qing-Lin Geng,? Min Luo," Yan Niu, An-Qing Wang,? Chen-Xu Wang,’ Meng 
Wang,’ Lei Zhang,” Liang Zhang,’ Rui-Kai Zhang,* Yao Zhang,* Ming-Gang Zhao,’ and Yang Zhou* 


! School of Physics and Microelectronics, Zhengzhou University, Zhengzhou 450001, China 
^ School of Physics, Nanjing University, Nanjing 210093, China 
*School of Physics, Nankai University, Tianjin 300071, China 
^ Institute of High Energy Physics, Chinese Academy of Sciences, 
19B Yuquan Road, Shijingshan District, Beijing 100049, China 
Research Center for Particle Science and Technology, 
Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao 266237, China 
S University of Chinese Academy of Sciences, 19A Yuquan Road, Shijingshan District, Beijing 100049, China 
School of Information Science & Engineering, Harbin Institute of Technology, Weihai 264209, China 
* School of Information Science and Engineering, Harbin Institute of Technology, Weihai 264209, China 


The reconstruction of the tracks of charged particles with high precision is crucial for HEP experiments to 
achieve their physics goals. The BESIII drift chamber, which is used as the tracking detector of the BESIII 
experiment, has suffered from aging effects resulting in degraded tracking performance after operation for ap- 
proximately 15 years. To preserve and enhance the tracking performance of BESIII, one of the proposals is to 
add one layer of a thin cylindrical CMOS pixel sensor based on state-of-the-art stitching technology between 
the beam pipe and the drift chamber. The improvement in the tracking performance of BESIII with such an 
additional pixel detector compared to that with only the existing drift chamber was studied using the modern 
common tracking software ACTS, which provides a set of detector-agnostic and highly performant tracking 
algorithms that have demonstrated promising performance for a few high-energy physics and nuclear physics 


experiments. 
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I. INTRODUCTION 


The Beijing Spectrometer (BESIII) [1] at the Beijing 
Electron-Positron (BEPCII) Collider has been a great success 
in producing promising physics results [2] in the 7-charm sec- 
tor since 2009. However, after operating for approximately 
15 years, the detectors at BESIII have been subjected to ag- 
ing effects, which degrade their performance [3]. 

The tracking detector at BESIII is a multilayer drift cham- 
ber (MDC), which provides measurements of the momentum 
and position of the charged tracks and information of energy 
loss in unit path length, i.e. dE/dx [4], of the charged tracks 
for particle identification. As shown in Ref. [3, 5], due to the 
beam-induced background with a hit rate up to 2 kHz/cm?, 
the gain of the MDC cells in the first ten layers has shown 
an obvious decrease, with a maximum decrease in approxi- 
mately 39% for the innermost layer cells in 2017. This further 
reduces the spatial and momentum resolution of the charged 
tracks. Because BESIII is not expected to complete its mis- 
sion in the foreseen years [2, 3], the track reconstruction per- 
formance must be preserved and enhanced to avoid compro- 
mising the physics goals of BESIII, along with a possible up- 
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grade of the tracking system with state-of-the-art technolo- 
gies for the detection of charged particles. To obtain good 
preparation for the potential malfunction of the MDC ow- 
ing to the aforementioned aging problem, plans to upgrade 
the BESIII inner tracker based on different technologies have 
been proposed [3]. This includes the replacement of the inner 
tracker with a new inner drift chamber [6] or a cylindrical gas 
electron multiplier (CGEM) tracker [7], which has attractive 
features such as high counting rate tolerance and low sensi- 
tivity to aging. 

Compared with gaseous detectors, silicon pixel detectors 
have excellent spatial resolution and good radiation resis- 
tance. Therefore, an additional option for the BESIII inner- 
tracker upgrade is to use a large-area thin complementary 
metal oxide semiconductor (CMOS) pixel sensor with good 
spatial resolution based on cutting-edge stitching technology, 
which has already been used to produce CMOS pixel sen- 
sors for medical imaging applications [8-10]. The design for 
a first wafer-scale stitched sensor prototype, known as the 
monolithic stitched sensor (MOSS) chip, is discussed in re- 
cent studies. These efforts aim to enhance the performance of 
the vertex detector, i.e. the ITS3, at the ALICE experiment, 
as detailed in Ref. [11]. 


A Common Tracking Software (ACTS) [12] is à common 
high-energy physics (HEP) software. It provides a set of 
detector-agnostic, high-performance, and modular tools for 
the track reconstruction in HEP. The software leverages mod- 
ern software technologies to facilitate concurrency, usabil- 
ity, maintenance, and extendability, addressing the prospec- 
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tive tracking challenges in HEP. To date, ACTS has been 
used for track reconstruction at ATLAS [13], FASER [14], 
sPHENIX [15], and STCF [16] for different types of tracking 
detectors. In particular, the promising tracking performance 
of ACTS for a tracking system with a drift chamber was first 
presented in Ref. [16]. 

In this study, the tracking performance of the BESIII MDC 
with an additional one-layer stitched cylindrical CMOS pixel 
detector inserted between the beam pipe and the inner wall of 
the MDC was studied using ACTS as the tracking software. 
The performance was compared with that of the current BE- 
SIII tracking detector using only the MDC based on the BE- 
SIII Offline Software System (BOSS) [17]. The remainder of 
this paper is organized as follows. In Section II, the BESIII 
MDC and the proposed pixel detector based on the cylindri- 
cal CMOS pixel sensor are presented. Section III introduces 
the tracking strategies in BOSS and ACTS. The improvement 
in the tracking performance at BESIII with an additional pixel 
detector is presented in Section IV. Finally, a brief conclusion 
is provided in Section V. 


Il. BESIII TRACKING DETECTOR AND STITCHED 
CMOS PIXEL DETECTOR 


The BESIII MDC is a cylindrical chamber operating with 
a helium-based gas mixture (He/C3Hg = 60:40) and is im- 
mersed in a 1-T magnetic field. The inner and outer radii of 
the MDC are 64 mm and 819 mm, respectively. The lengths 
of the wires range from 774 mm for the innermost layer to 
2400 mm for the outermost layer. The drift cells, which are 
almost square-shaped, are arranged in 43 circular layers and 
alternate between the stereo layers and axial layers, i.e. in the 
order of eight stereo layers, 12 axial layers, 16 stereo layers, 
and 7 axial layers. The cell dimensions are approximately 
12 mm x 12 mm for the eight inner layers and 16.2 mm x 
16.2 mm for the 35 outer layers. The inner chamber is com- 
posed of 8 inner layers, and the outer chamber is composed 
of 35 outer layers. 

CMOS pixel sensors have been widely used for vertex de- 
tectors in HEP because of their excellent spatial resolution 
down to a few um, tolerance of high hit rates up to 10° 
Hz/cm?, and good detection efficiency and radiation resis- 
tance. For example, the CMOS pixel sensor has been used 
as the inner tracker or vertex detector at STAR experiments 
[18], ALICE experiments [19] and sPHENIX experiments 
[20]. It is also considered as the primary choice for pixel 
sensor technology in the CEPC physics program [21]. De- 
spite the advantages of the CMOS pixel sensor, the tracking 
resolution of a traditional CMOS pixel detector, particularly 
for low-momentum tracks, can be limited by the irreducible 
material budget of the support structure and cooling pipe. 
This was demonstrated in Ref. [22], in which a MAPS-based 
CMOS pixel sensor was proposed to upgrade the BESIII in- 
ner chamber. Therefore, tracking based on traditional CMOS 
pixel sensors can be challenging at BESIIL in which tracks 
with a transverse momentum as low as 150 MeV must be re- 
constructed with good resolution. 
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In recent years, innovative stitching technologies have 
emerged with the aim of producing large-area sensors on one 
wafer, which can be further thinned to approximately 50 jum. 
With this thickness, the wafer can be curved into a cylindrical 
shape, thereby vastly simplifying the support structure and 
reducing the material budget. Based on this stitching technol- 
Ogy, it is possible to construct a vertex detector using large- 
area thin CMOS pixel sensors with cylindrical shapes and 
simplified support structures. As discussed in Ref. [11], the 
material budget per pixel layer is expected to be reduced to 
0.05% Xo to allow bending of the silicon. In this study, we 
considered a vertex detector comprising a single layer of a 
stitched cylindrical pixel sensor with a thickness of 50 wm 
(0.05% Xo), as proposed in Ref. [11], length of 220 mm (in- 
cluding an insensitive area within 7 mm of each end), and 
a resolution of 8.66 um and 57.74 um in the r-ó and z di- 
rections, respectively. This layer was inserted between the 
63-mm diameter beam pipe [1] and the MDC inner wall. The 
stitching technology foresees cooling by airflow and no read- 
out circuits in the active area. Therefore, the additional ma- 
terial budget from the cooling and readout infrastructure, as 
well as the support structure in the active region of the pixel 
layer, is expected to have a negligible impact on the tracking 
performance; hence, it was not considered in this study. 

The geometry of the BESIII MDC and the pixel detector is 
shown in Figure 1. 


Fig. 1. Clipped view of the geometry of the BESIII MDC (cell wires 
are shown in green colors) and the pixel detector (shown in orange 
color) outside of the beampipe. 


Ill. TRACKING WITH BESIII OFFLINE SOFTWARE 
AND ACTS 


The simulation study was based on samples generated us- 
ing BOSS, an offline software framework for event generation, 
simulation, reconstruction, performance validation, physics 
analysis, and visualization [23] for the BESIII experiment. 

In BOSS, the production of both charmonium resonances 
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and the continuum process from et e~ collisions is provided 
by the KKMC [24] generator, and the decays of particles 
were modeled using EVTGEN [25], which can also be used 
to generate single-particle samples. The detector descrip- 
tion is based on the Geometry Description Markup Language 
(GDML) [26], and the interaction of the particles with the de- 
tectors is simulated using GEANTA [27]. 

The performance of the BESIII MDC with the additional 
pixel detector was studied using the ACTS software and com- 
pared with the performance of the BESIII MDC studied using 
the BESIII tracking software within BOSS. The workflow for 
studies of the tracking performance is illustrated in Figure 2. 


Event generation Magnetic 
field map 


ACTS tracking 
geometry 


MDC+Pixel 
measurements 
in ACTS 


Í MDC +Pixel / 
; geometry : 


Track finding 


ACTS tracking 
with CKF 


Tracks 


Fig. 2. Workflow of tracking performance studies based on BOSS 
and ACTS. 


A. Tracking in BOSS 


The trajectory of a charged particle in a magnetic field is 
parameterized using the helix track parameters at BESIII, as 
described in Ref. [28]. Track finding and track fitting are 
two tasks in BOSS track reconstruction. Track finding is a 
pattern-recognition problem of classifying measurements into 
subsets and creating track candidates, whereas track fitting 
is an estimation of the helix parameters. Two basic track- 
finding algorithms were implemented in BOSS: the template- 
matching algorithm (PAT) [29] and the track segment-finder 
algorithm (TSF) [28]. A specialized track-finding method 
called TCurlFinder [30] and a global track-finding method 
based on the Hough transform (HOUGH) [31] have been im- 
plemented to salvage low transverse momentum tracks with 
pr < 120 MeV. In addition, an extended segment con- 
struction scheme that can achieve higher efficiency for low- 
transverse-momentum tracks has been developed [32]. The 
tracks found using these algorithms are combined and fed into 
a Kalman filter [33] for track fitting. Further track extrapola- 
tion is performed to obtain the track parameters of the other 
subdetectors at BESIII. 
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B. Tracking with ACTS 


A detailed introduction to ACTS, including the geom- 
etry description, parameterization of the track parameters 
and measurements and tracking algorithms, can be found in 
Ref. [12]. The implementation of ACTS for track reconstruc- 
tion of the BESIII tracking detector is outlined below. 

To use ACTS for track reconstruction, the detector geom- 
etry with a detailed description of the detector material and 
placement, e.g. the GEANT4-based detector geometry, must 
be transformed into ACTS internal geometry, i.e. ACTS track- 
ing geometry. This geometry has a simplified description of 
the passive detector material to facilitate fast navigation and 
track reconstruction. The TGEO [34] plugin in ACTS was used 
to transform the TGEO version of the BESIII detector geom- 
etry, which was created based on GDML exported from the 
GEANT4-based geometry, into the ACTS tracking geometry 
with material mapped to the auxiliary surfaces. The one pixel 
layer of the pixel detector was converted to a cylinder surface 
in ACTS and the 43 layers of signal wires of the MDC were 
converted to 43 layers of line surfaces, dedicated to represen- 
tation of measurements for a drift chamber or track parame- 
ters at the interaction point? in ACTS. 

The hits on the MDC after digitization were transformed 
into one-dimensional measurements associated with the line 
surfaces in ACTS, describing the drift distance of the drift 
chamber. The hits on the pixel detector after the simula- 
tion were transformed into two-dimensional measurements 
by smearing the simulated hits with the resolution of the pixel 
detector using Gaussian functions. These measurements are 
associated with ACTS cylinder surfaces, in which the local x 
(y) coordinate represents r - y (z) in the cylinder frame; that 
is, r represents the radius of the cylinder, y represents the 
azimuthal angle of the position on the cylinder, and z is the 
coordinate in the z direction. 

The BESIII magnetic field was transformed into an ACTS- 
interpolated magnetic field using the field map of the BESIII 
magnetic field. With an interpolated field provider in ACTS, 
the value of the magnetic field for any given position is calcu- 
lated by interpolating from a grid of known values, e.g. eight 
corner points of a field cell in three-dimensional coordinate 
system. 

Track fitting and finding are performed simultaneously us- 
ing the combinatorial Kalman filter [35] algorithm in ACTS, 
which has the capability of rejecting noise hits based on the 
X? calculated using the distance between the hit and the pre- 
dicted track parameters and the covariance of the hit and the 
predicted track parameters. If multiple hits are found to be 
compatible with the predicted track parameters, the hit with 
the best x? is used to filter out the track parameters. The 
smoothed track parameters from the first measurement are 
extrapolated to the beam line to obtain the estimated track 
parameters at the interaction point. 


3 Suppose the wire has a direction of w and the track direction is ¢ in the 
global coordinate frame, the x axis and y axis in the local coordinate frame 
of a line surface is W x t and Ù, respectively. 
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IV. TRACKING PERFORMANCE STUDIES 
A. Monte-Carlo sample generation 


The effect of the pixel detector on the resolution of the 
track parameters was studied using single u~ and single m~ 
samples. The samples were generated with a fixed trans- 
verse momentum pr, cos (0 is the polar angle) uniformly 
distributed between [-0.8, 0.8], and an azimuthal angle $ uni- 
formly distributed in the range of [0, 27]. The impact of the 
pixel detector on the tracking efficiency was studied using 
(3686) — «m Jf, J/ — pt pm events generated at 
y's = 3.686 GeV. The two-dimensional distributions of cosó 
and pr for u and 7 during this process are shown in Figure 3. 

Random background hits in the MDC due to beam-related 
background or electronic noise were generated in the stan- 
dard simulation sample production at BESIII. Because a ded- 
icated background noise model for the pixel detector is not 
yet available, no background noise was considered for the in- 
serted pixel layer in this study. 


B. Track reconstruction performance 


The resolution of the impact track parameters, do and zo, 
and the relative resolution of the transverse momentum pr 
as a function of pr for the BESIII MDC with an additional 
pixel layer, which is placed with radius rpjxe) at three different 
values (i.e. 35 mm, 45 mm, or 55 mm) compared to those of 
BESIII MDC, are shown in Figure 4. Although rpjxe) has little 
impact on the resolution of pr, the spatial resolution is better 
with a smaller rpixe1. With rye; = 35 mm, the resolution of 
do and zo is 60 um and 120 um for u` and 7^ with pp = 1 
GeV, respectively, and the relative resolution of 44% for pr 
at pr = 1 GeV was achieved. The resolutions of do, zo, pr 
can be improved by up to 67%, 93%, and 32%, respectively, 
by adding the proposed pixel detector to the current BESIII 
tracking detector. 

A comparison of the tracking efficiency, which is defined 
as the fraction of the generated charged particles which have 
the corresponding reconstructed tracks in the detector region 
(Icos0| < 0.93), between the BESIII MDC with the addi- 
tional pixel layer and the BESIII MDC alone is shown in Fig- 
ure 5. The additional pixel layer can significantly improve 
the tracking efficiency for tracks with a particularly low mo- 
mentum or large Icos6l. For u with pr below 0.6 GeV in 
the (3686) > ntr J/p, J/ — wtp process, the im- 
provement can reach 996. The improvement can reach 1096 
for 7 with 0.85 < Icos6| < 0.93. A pixel with ry, = 35 mm 


[1] M. Ablikim, et al, Design and construction of the BE- 
SIII detector, Nucl. Instrum. Meth. A 614 (3) (2010) 
345-399. doi:https://doi.org/10.1016/j.nima. 
2009.12.090. 
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provides better tracking efficiency than the other two choices. 


V. CONCLUSION 


After operating for approximately 15 years, the tracking 
detector of BESIII experiment has suffered from aging ef- 
fects and must be upgraded to preserve its tracking perfor- 
mance for BESIII to fullfill its remaining physics goals in the 
next few years. The possibility of inserting a one-layer pixel 
detector using a large-area thin cylindrical CMOS pixel sen- 
sor based on stitching technology between the beam pipe and 
the MDC was considered. The spatial and momentum resolu- 
tions were studied using ACTS and compared with the perfor- 
mance of the current BESIII tracking detector obtained using 
the BESIII offline software. The proposed method based on 
the stitched cylindrical CMOS pixel sensor was found to be 
promising, i.e. the resolution of do, zo, pr can be improved 
by up to 67%, 93%, and 32%, respectively, and the tracking 
efficiency for particles in the (3686) — m ^m J/w, J/v > 
L^ events can be improved by up to 10% in the low pr 
or large lcos0l region, with the additional pixel layer placed 
at a radius of 35 mm. Meanwhile, this is the first time that 
ACTS has been used for track reconstruction for a drift cham- 
ber with calibrated measurements and realistic noise hits. In 
future studies, the noise model of this state-of-the-art pixel 
detector will be studied, and the tracking performance con- 
sidering the noise of the pixel detector will be investigated. 
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Fig. 5. The tracking efficiency of ji (top), m (bottom) in the process of 4»(3686) — ata J/w, J/ — wtp as a function of particle pr 
(left) and cos0 (right) for the BESIN MDC only (black circle), and BESIII MDC with an additional pixel layer (denoted as "Pixel") placed 
with rpixe: = 35 mm (blue dot), 45 mm (purple triangle) and 55 mm (yellow triangle), respectively. 
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